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We experimentally demonstrate the coherent oscillations of a tunable superconducting flux qubit 
by manipulating its energy potential with a nanosecond-long pulse of magnetic flux. The occupation 
probabilities of two persistent current states oscillate at a frequency ranging from 6 GHz to 21 GHz, 
tunable via the amplitude of the flux pulse. The demonstrated operation mode allows to realize 
quantum gates which take less than 100 ps time and are thus much faster compared to other 
superconducting qubits. An other advantage of this type of qubit is its insensitivity to both thermal 
and magnetic field fluctuations. 

PACS numbers: 03.67.Lx, 85.25.Dq 
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Superconducting qubits stand between the most 
promising systems for the realization of quantum com- 
putation. Coherent quantum evolution and manipula- 
tion hav_e_been_de_monstrated and extensively studied for 
coupled superconducting qubits 
most cases, the state of supercon- 
ducting qubits are manipulated by means of microwave 
pulses, with a technique similar to the NMR manipula- 
tion of atoms. An alternative way to manipulate qubits 
is based on modifying their energy potential without ap- 
plying any microwave signals . The latter approach 
requires a much simpler experimental technique and of- 
fers the possibility of using classical logic signals to con- 
trol a quantum processor in situ, which is advantageous 
for the large scale implementation of a quantum circuits. 

In this Letter, we report the observation of tunable co- 
herent oscillations in a SQUID-based flux qubit. These 
oscillations are obtained by manipulating the qubit with 
nanosecond-long pulses of magnetic flux rather than mi- 
crowaves. By this technique, we could increase the oscil- 
lation frequency up to 21 GHz, which allows to perform 
very fast logical quantum gates. Since the relevant qual- 
ity factor of a qubit is the number of gate operations 
which can be performed during its coherence time, this 
result is of particular interest towards the realization of 
a solid-state quantum computer. 

The investigated circuit, shown in Fig. [Ija), is a dou- 
ble SQUID consisting of a superconducting loop of induc- 
tance L — %b pH, interrupted by a small dc SQUID of 
loop inductance / — 6 pH. This dc SQUID is operated as 
a single Josephson junction (JJ) whose critical current 
is tunable by an external magnetic field. Each of the 
two JJs embedded in the dc SQUID has a critical cur- 
rent /o — 8^A and capacitance C — 0.4 pF. The qubit is 
manipulated by changing two magnetic fluxes and 
applied to the large and small loops by means of two coils 



of mutual inductance = 2.6 pH and Mc = 6.3 pH, 
respectively. The readout of the qubit flux is performed 
by measuring the switching current of an unshunted dc 
SQUID, which is inductively coupled to the qubit [l^ . 
The circuit was manufactured by Hypres [ist using stan- 
dard Nb/AlOx/Nb technology in a 100 A/cm^ critical 
current density process. The dielectric material used for 
junction isolation is Si02. The whole circuit is designed 
gradiometrically in order to reduce magnetic flux pick-up 
and spurious flux couplings between the loops. The JJs 
have dimensions of 3 x 3/im^ and the entire device occu- 
pied a space of 230 x 430/im^. All the measurements have 
been performed at a sample temperature of 15 mK. The 
currents generating the two fluxes and $c were sup- 
plied via coaxial cables including 10 dB attenuators at 
the IK-pot stage of a dilution refrigerator. To generate 
the flux a bias-tee at room temperature was used to 
combine the outputs of a current source and a pulse gen- 
erator. For biasing and sensing the readout dc SQUID, 
we used superconducting wires and metal powder filters 
14| at the base temperature, as well as attenuators and 



low-pass filters with a cut-off frequency of 10 kHz at the 
IK-pot stage. The chip holder with the powder filters 
was surrounded by one superconducting and two cryop- 
erm shields. 

Assuming identical junctions and negligible inductance 
of the smaller loop {I ^ L), the system dynamics is 
equivalent to the motion of a particle with the Hamil- 
tonian 
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where (p = ^/^b is the spatial coordinate of the equiv- 
alent particle, p is the relative conjugate momentum, 
M — C^l is the effective mass, (fx = ^x/^b and 
tpc — are the normalized flux controls, and 
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FIG. 1: (a) Schematic of the flux qubit circuit, (b) The 
control flux "I>c changes the potential barrier between the two 
flux states \L) and |_R), here = 0.5 $0. (c) Efltect of the 
control flux on the potential symmetry. 



P{(Pc) = {2IoL/^h) COS ipc, with $0 = h/ {2e) and $fc 
$0/ (27r). For /3 < 1 the potential has a single minimum, 
otherwise it consists of multiple wells. In the particular 
case of 1 < /3 < 4.6 and = $o/2, the system potential 
is a symmetric double well shown in Fig. [Ijb). The two 
states \L) and \R), which are respectively localized in the 
left and right potential well, correspond to a persistent 
current circulating either clockwise or counter-clockwise 
in the main SQUID loop. As it is shown in Fig.[Tl^b), the 
external flux <I>c controls the height of the barrier sepa- 
rating the minima, while a variation of changes the 
symmetry of the potential as indicated in Fig. [Ijc) . In 
this work, we exploit both the double well and the single 
well properties. The double well potential shape is used 
for qubit initialization and readout. The single well, or 
more exactly the two lowest energy states |0) and |1) in 
this well, is used for the coherent evolution of the qubit. 

We use a well established procedure [l^ to identify 
the regions where the system has a double well poten- 
tial in the $c - plane. The flux response $ of the 
qubit is measured as a function of and $c fluxes 
and the switching points between different flux states 
are detected. Figure ^a,) shows $ - $3, characteristics 
obtained for two values using initial state prepara- 
tion in different wells. One can easily identify here a 
bi-stability region in the vicinity of ^x ~ 0.5 $0 and a 
mono-stability region outside the hysteretic curve, cor- 
responding to the presence of a double and a single po- 
tential well, respectively. At the border between these 
regions we find abrupt switching between the two sta- 
ble states. The positions ^x of the switching points are 
plotted in Fig.[2{b) for different $c values. This diagram 
allows to easily identify the combinations of parameters 
resulting in a single- or double-well potential. A single- 
well region is found for $c ^ 0.42. 

The measurement process that we used to observe 
coherent oscillations consists of several steps shown in 
Fig. [2Ia). Each step is realized by applying a combina- 
tion of magnetic fluxes ^x and <I>c as indicated by num- 
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FIG. 2: (a) The measured double SQUID flux ^ in depen- 
dence of plotted for two different values of $c and ini- 
tial preparation in either potential well, (b) Position of the 
switching points (dots) in the <l?c — $a; parameter space. Num- 
bered tags indicate the working points for qubit manipulation 
at which the qubit potential has a shape as indicated in the 
insets. 



bers in Fig. (H^b). The first step in our measurement is 
the initialization of the system in a defined flux state 
(1). Starting from a double well at ^x — ^o/2 with 
high barrier, the potential is tilted by changing ^x until 
it has only a single minimum (left or right, depending 
on the amplitude and polarity of the applied flux pulse). 
This potential shape is maintained long enough to ensure 
the relaxation to the ground state. Afterwards the po- 
tential is tuned back to the initial double- well state (2). 
The high barrier prevents any tunneling and the qubit 
is thus initialized in the chosen potential well. Next, 
the barrier height is lowered to an intermediate level (3) 
that preserves the initial state and allows to use just a 
small-amplitude $c flux pulse for the subsequent manip- 
ulation. The following $c-pulse transforms the potential 
into a single well (4). The $c-pulse duration At is in the 
nanosecond range. The relative phase of the ground and 



the first excited states evolves depending on the energy 
difference between them. Once $c-pulse is over, the dou- 
ble well is restored and the system is measured in the 
basis {\L), \R)} (5). The readout of the qubit flux state 
is done by applying a bias current ramp to the dc SQUID 
and recording its switching current to the voltage state. 

The pattern that realizes the above described manip- 
ulation is reported in Fig. [3Jb) . The flux is switched 
between two values: is used to create a strongly 

asymmetric potential for qubit initialization in the left 
or right well, and equal to $o/2 (or very close to 
it) transforms the potential into a symmetric (or nearly 
symmetric) double well. The flux <I>c is changed between 
three different values: $ci and $c2 define, respectively, 
high and intermediate amplitudes of the barrier between 
the two minima, while at = ^c2 + A$c the barrier is 
removed completely and the potential turns into a single 
well. The amplitude A$c of the pulse is varied allowing 
for single wells of different curvature at the bottom. The 
nominal rise and fall times of this pulse are t^/ / = 0.6 ns. 

The flux pattern is repeated for 10^ — 10'' times in 
order to evaluate the probability Pl = |(i|^finai)|^ of 
occupation of the left state at the end of the manipula- 
tion. By changing the duration At of the manipulation 
pulse <i>c, we observed coherent oscillations between the 
occupations of the states \L) and \R) shown in Fig. SJa). 
The oscillation frequency could be tuned between 6 and 
21 GHz by changing the pulse amplitude A<I>c. These 
oscillations persist when the potential is made slightly 
asymmetric by varying the value ^xi- As it is shown in 
Fig-IHb), detuning from the symmetric potential by up to 
±2.9m$o only slightly changes the amplitude and sym- 
metry of the oscillations. When the qubit was initially 
prepared in |_R) state instead of \L) state we observed 
similar oscillations. 

To understand the physical process behind the ob- 
served oscillations, let us discuss in detail what hap- 
pens during the manipulation. Suppose the system is 
initially prepared in the left state \L) of a perfectly sym- 
metric double well potential. During the $c pulse, the 
potential has only one central minimum and can be ap- 
proximated by a harmonic oscillator potential with fre- 
quency cJo(*c3) « l/V2lC^l-/3($c3). The pulse 
transforms the initially prepared left state (that is a 
symmetric superposition of the two lowest energy eigen- 
states of the double well potential |0) and |l), i.e. \L) = 
(|0) -I- |i))/V2 ) into the superposition of the two lowest 
energy eigenstates |0) and |1) of the single- well poten- 
tial. To achieve that, the pulse rise time needs to be 
shorter than the relaxation time but, at the same time, 
long enough to avoid population of upper energy lev- 
els. During the plateau of the A<I>c pulse, the relative 
phase 9 between the states |0) and |1) evolves with time 
at the Larmor frequency given by too = (-^'i — Eq) /h. 
At the end of the pulse the accumulated relative phase 
becomes 9 = wpAt. Turning the system back into the 
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FIG. 3: (a) Variation of the potential shape during the manip- 
ulation, (b) Time sequence of the readout dc SQUID current 
(topmost line) and flux bias values (bottom lines). 



double- well maps the phase to the two flux states \L) 
and \R). The final state after the flux pulse <I>c is 
|*finai) = cos(6l/2) \L) +ism{9/2) \R). Note that in the 
more realistic case of non perfectly symmetric double- 
well potential, the initial left state is no more a symmet- 
ric superposition, but tends more to either |()) or |i^ due 
to the potential unbalancing. However, a pulse with a 
short rise time induces a non-adiabatic transition that 
populates mainly the two lowest energy eigenstates |0) 
and |1) in the single well potential. This condition can 
be met in a narrow region of the flux bias plane called 
"portal" [l^ . This non-adiabatic transition also leads to 
the phase evolution process described above. 

In order to verify above interpretation, we numerically 
solved the time-dependent Schrodinger equation for this 
system. The simulation showed that with our experimen- 
tal parameters the transition between the first two levels 
occurs as described, while the occupation of upper levels 
remains below few percents. 

The oscillation frequency luq depends on the ampli- 
tude of the manipulation pulse since it determines 
the shape of the single well potential and the energy 
level spacing Ei — Eq. A pulse of larger amplitude A<i>c 
generates a deeper well having a larger level spacing, 
which leads to a larger oscillation frequency as shown 
in Fig. mja). In Fig. El we plot the energy spacing be- 
tween the ground state and the three excited states (in- 
dicated as {Ek — Eq) /h with k=l,2,3) versus the flux 
<I>c3 ~ $c2 + A$c obtained from a numerical simulation 
of our system using the experimental parameters. In the 
same flgure, we plot the measured oscillation frequencies 
for different values of <I>c (open circles). Excellent agree- 
ment between simulation (solid line) and data strongly 
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FIG. 4: Probability to measure the state \L) in dependence 
of the pulse duration At for the qubit initially prepared in 
the \L) state, and for (a) different pulse amplitudes A$c, 
resulting in the indicated oscillation frequency, and (b) for 
different potential symmetry by detuning $x from $o/2 by 
the indicated amount. 



supports our interpretation. The fact that a small asym- 
metry in the potential does not change the oscillation 
frequency, as shown in Fig.[4ljb), is consistent with the in- 
terpretation as the energy spacing Ei — Eq is only weakly 
affected by small variations of This provides protec- 
tion against noise in the controlling flux $2,. 

The measured oscillation decay time of about 2 ns is of 
the same order as the Rabi oscillations decay and the en- 
ergy relaxation time Ti which we measured by using stan- 
dard microwave 7r-pulse manipulation on the same de- 
vice. It is also comparable to the coherence time obtained 
on similar devices fabricated using the same technology 
17| , suggesting that coherence is not limited by the ma- 
nipulation procedure reported in this paper. We believe 
that the decay time of the reported high-frequency coher- 
ent oscillations can be increased by two orders of magni- 
tude by reducing the area of the JJs and using an appro- 
priate dielectric instead of Si02 as insulating material in 
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FIG. 5: Calculated energy spacing of the first (solid line), sec- 
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respect to the ground state in the single well potential, plotted 
vs. the control flux amplitude $c3. Circles are the experimen- 
tally observed oscillation frequencies for the corresponding 
pulse amplitudes. 



the junction fabrication [18| . 

We note that a similar qubit manipulation procedure 
has been reported by Koch et al. , demonstrating Lar- 
mor oscillations in a flux qubit coupled to a harmonic 
oscillator. It should be emphasized that, in our case, the 
oscillator is not required, which simplifies the realization 
of the qubit circuit. Moreover, in contrast to Ref. 0], 
our approach provides a wide range tunability of the fre- 
quency of coherent oscillations. This allows for a variety 
of quantum gates by manipulating both duration and 
amplitude of the flux pulses, and moreover is important 
concerning the realization of a controllable coupling be- 
tween qubits or quantum busses such as resonant cavities. 

In conclusion, we presented the coherent manipulation 
of a flux qubit without using microwaves. The reported 
approach seems particularly promising for the realization 
of circuits with many qubits, and it appears to be well 
suited for integration with RSFQ control electronics [l^ . 
The benefits of the reported system are the possibility 
of in situ tuning the frequency of oscillations and their 
insensitivity to small changes in the potential symmetry. 
The high frequency of oscillations allows for very fast 
qubit gate operations, and the large energy gap between 
the qubit states during coherent evolution protects the 
system from thermal activation to upper energy states. 
Moreover, the oscillation frequency depends only weakly 
on the control pulse amplitude A$cj in contrast to the 
exponential sensitivity of the oscillations in a double well 
potential [20] , which makes the qubit manipulation more 
reliable. 
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